Draft version July 19, 2012 

Preprint typeset using I^T^jK style emulateapj v. 08/22/09 



REDSHIFTS, SAMPLE PURITY, AND BCG POSITIONS FOR THE GALAXY CLUSTER CATALOG FROM 
THE FIRST 720 SQUARE DEGREES OF THE SOUTH POLE TELESCOPE SURVEY 

J. Song\ a. Zenteno^'', B. Stalder", S. Desai^ '\ L. E. Bleem''-''. K. A. Aird'', R. Armstrong*, M. L. N. Ashby*'^, 
M. Bayliss'', G. Bazin^-\ B. A. Benson"'^", E. Ber.tin^\ M. Brodwin", J. E. Carlstrom^''''1°'1=''", C. L. Ciiang^'^"'", 
H. M. Cno^\ A. Clocchiatti"', T. M. Crawford^-^-', A. T. Crites'"''^'', T. de Haan^^ M. A. Dobbs^', J. P. Dudley", 
R. J. FoleyS E. M. GEORGE^^ D. Gettings"', M. D. Gladders"-^\ a. H. Gonzalez"', N. W. Halverson^", 
N. L. Harrington^*, F. W. Higii^^'\ G. P. Holder", W. L. HoLZAPFEL^^ S. Hoover-'"', J. D. HRUBES^ M. Joy^^, 
R. KEISLER=■^ L. Knox^^ a. T. LEE^*'2^ E. M. Leitch^'", J. LIU^'^ M. Lueker^^'^*, D. Luong-Van^ 

D. P. MARRONE2^ M. McDONALD^^ J. J. McMAHOn\ J. MBHL^'^^ , S. S. MEYER^''''!"'", L. MOCANU^'l*, J. J. Mohr2'3 ^7^ 

T. E. MONTROY^*, T. NATOLI^'^, D. NURGALIEv'', S. PADIN^'^^'^* T. PLAGGE^'^^, C. PRYKE^^ C. L. Reichardt^* 
a. REST^°, J. RUEL^ J. E. RUHL^^ B. R. SaLIWANCHIK^* , a. SaRO , J. T. SAYRE^^ K. K. SCHAFFER^'1'''^\ L. SHAW"•^^ 

E. Shirokoff^*^-^", R. Suiiada^ H. G. Spieler^'', S. A. Stanford^^-^-'' , Z. Staniszewski^*^, A. A. Stark", K. Story^'^, 
C. W. Stubbs"'", a. van Engelen", K. Vanderlinde", J. D. Vieira'-*''^*, R. Williamson'"'"', and O. Zahn^'''^'* 

Draft version July 19, 2012 

ABSTRACT 

Wc present the results of the ground- and space-based optical and near-infrared (NIR) follow-up of 
224 galaxy cluster candidates detected with the Sunyaev-Zel'dovich (SZ) effect in the 720 deg^ of the 
South Pole Telescope (SPT) survey completed in the 2008 and 2009 observing seasons. We use the 
optical/NIR data to establish whether each candidate is associated with an overdensity of galaxies 
and to estimate the cluster redshift. Most photometric redshifts are derived through a combination 
of three different cluster redshift estimators using red-sequence galaxies, resulting in an accuracy of 
Az/(1 + z) = 0.017, determined through comparison with a subsample of 57 clusters for which we 
have spectroscopic redshifts. We successfully measure redshifts for 158 systems and present redshift 
lower limits for the remaining candidates. The redshift distribution of the confirmed clusters extends 
to 2; = 1.35 with a median of z^^.^ = 0.57. Approximately 18% of the sample with measured redshifts 
lies at 2; > 0.8. We estimate a lower limit to the purity of this SPT SZ-selected sample by assuming 
that all unconfirmed clusters are noise fluctuations in the SPT data. We show that the cumulative 
purity at detection significance ^ > 5 ( ^ > 4.5) is > 95% (> 70%). Wc present the red brightest 
cluster galaxy (rBCG) positions for the sample and examine the offsets between the SPT candidate 
position and the rBCG. The radial distribution of offsets is similar to that seen in X-ray-selected 
cluster samples, providing no evidence that SZ-selected cluster samples include a different fraction of 
recent mergers than X-ray-selected cluster samples. 

Subject headings: galaxies: clusters: general — galaxies: distances and redshifts — cosmology: obser- 
vations 
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1. INTRODUCTION 



In November 2011, the South Pole Telescope (SPT; 



Carlstrom et al. 2011 1 collaboration completed a 2500 
deg survey, primarily aimed at detecting distant, mas- 
sive galaxy clusters thro ugh their Sunyae v -Zel'd ovich 
(SZ) effect signature. In [Reichardt et al.| ( |2012[ R12 
hereafter) , the SPT team presented a catalog ot 224 clus- 
ter candidates from 720 deg^ observed in the 2008-2009 
seasons. In this work, we present the optical and near- 
infrared (NIR) follow-up observations of the cluster can- 
didates reported in R12, mainly focusing on follow-up 
strategy, confirmation and empirical purity estimate for 
the cluster candidates, photometric redshift estimations 
of confirmed clusters, and the spatial position of the red 
brightest cluster galaxies. 

Galaxy clusters have long been used f or the study of 



et al. 



Geller & Beers 



1982 


White et al. 


1993 


I. Soon after the discovery ot 


the cosmic acceleration ( 


Schmidt et al. 1998 


Perlmutter 



redshitt evolution of the cluster mass function could pro 
vide a p owerful tool to further und erstan d the underlying 
causes (Wang fc Steinhardt||1998 Haiman et al. 2001 



Holder eriTp OOl ; Battyc Wcller 2001 

cise the oretical investigations (Majumdar 

Hul[200 3; Majumdar & Mohr 2004; Molnar eTaT] 



More pre- 



Mohr 



2003 



Wang et al.||2004[ |Lima fc Hu|2005( |2007| ) identified the 
key challenges associated with cluster surveys, which in- 
clude: (1) producing large uncontaminated samples se- 
lected by an observable property that is closely related to 
the cluster mass, (2) measuring cluster redshifts for large 
samples and (3) precisely calibrating the cluster masses. 

Competitive approaches to producing large cluster 
samples include optical multiband su rveys (e.g. 



ders fc Yeel|2005 } [Koest er et al. 
'tiscnhardt et al. 



povichi 
Pacaui 



, X-ray s urveys 



1 



Glad- 



2007), infrar ed surveys 
uzzin et ai.l 2008; PaP 



e-g 



I'inoguenov et al. 2007 



eTaLllMTl |Vikhlinin et al.||:^0(jyi |i-'inoguenov 



et al.||20101 |Mantz et al.||20101 [Fassbender et al.|120rii 



millimeter- wave 
2010, 



Lloyd-Davies et al. 2011 Suhada et ah 2012), and 



mm- wave) surveys ( Vanderlindc et al 



Planck Collaboration et al 



Marriag e et al.[|2011 
2011| [Reichard t et al.|2012p. The mm- wave surveys cap 
italize on the cluster SZ effect signature, which is pro 



duced by the inverse Compton scattering of cosmic mi- 
crowave backgrou nd photons by the energet ic electrons 
within the cluster ( Sunyaev fc Zerdovichh972| ) . The sur- 
face brightness of the SZ effect is redshitt-independent, 
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making SZ surveys a particularly powerful tool for iden- 
tifying the most distant clusters. It is typical for X-ray 
and mm-wave surveys to have accompanying multiband 
optical imaging to enable photometric redshift measure- 
ments; these multiband optical data also enable a second 
stage of cluster candidate confirmation, verifying the pu- 
rity estimation of the X-ray or SZ-selected cluster sam- 
ples. 

Ideally, one would coordinate an SZ survey with a deep, 
multiband optical survey over the same region; indeed, 
the Da rk Energy Survey (DES^^ Cease et al. 2008^, Mohr] 
et al.(20 08) and the SPT are coordinated m tins way. Be- 
cause of the different development timelines for the two 
projects, it has been necessary to undertake extensive 
cluster-by-cluster imaging follow-up for SPT using a se- 
ries of ground-based telescopes together with space-based 
NIR imaging (from Spitzer and WISE) . The NIR data are 
of particular importance in the confirmation and redshift 
estimation of the z > 1 massive galaxy clusters, which are 
especially interesting for both cosmological studies and 
studies of the evolution of clusters themselve s. Pointed 
observations were used in High et al. ( [2010 HIO here- 
after) to pro vide redshift and richness estimates of the SZ 
detections of Vandcrhnde et al. ( 2010^, and subsequ ently 
by iWilliamson et al., ( .20TT| and ptory et al.| ( |2011 1 . 

Cluster samples from high-resolution SZ surveys can 
be also used to explore the evolution of cluster prop- 
erties as a function of redshift. Previous studies using 
X-ray-selected clusters have identified a correlation be- 
tween the dynamical state of a cluster and the projected 
offset between t he X-ray centers and the brightest clus- 
ter galaxy (e.g., Katayama et al. '120031, 'Sanderson_et al. 



2009 



ga laxy (e.g., [Katayama et a,l."Z[)[)6\ 'banderson etaL 
9] Haarsma et al. 2010| Eassbendcr et al. 20lTj). In 



princi]3le, SZ-selected clusters can serve as laboratories to 
search for this correlation also, if the spatial resolution of 
SZ detections is high enough to detect the significance of 
offsets between the SZ centers and the brightest cluster 
galaxies. Systematic comparison between X-ray and SZ 
samples will indicate if the selection of the two methods 
differs in terms of the dynamical state of clusters. 

This paper is structured as follows: we briefly describe 
the SPT data and me thods for extracting the cluster 
sample in |2.1[ In j ]2.2[ we provide details of the follow- 
up strategy, as well as data processing. ^ is dedicated 
to a detailed description of the analysis of our follow- 
up data, including redshift estimation using optical and 
Spitzer data, the derivation of redshift lower limits for 
those systems that are not confirmed, and the selection 
of red brightest cluster galaxies (rBCGs) in the clusters. 
Results are presented in ^ and discussed further in SjSj 
Throughout this paper, we use the AB magnitude system 
for optical and NIR observations unless otherwise noted 
in the text. 



2. DISCOVERY & FOLLOWUP 

2.1. SPT Data 

Here we briefly summarize the analysis of the SPT 
data and the extraction of cluster candidates from that 
data; we refer th e reader to R12 and previous SPT clus- 
ter pu blica tions ([Staniszewski et al 2009 Vanderlinde 
et al.||2010| iWilhamson et al.||2011[ ) Tor more details. 



www . darkenergy survey, org 
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The SPT operates in three frequency bands, although 
only data from the 95 GHz and 150 GHz detectors were 
used in finding clusters. The data from all detectors at 
a given observing frequency during an observing period 
(usually 1-2 hours) are combined into a single map. The 
data undergo quality cuts and are high-pass filtered and 
inverse-noise- weighted before being combined into a map. 
Many individual-observation maps of every field are co- 
added (again using inverse noise weighting) into a full- 
depth map of that field, and the individual-observation 
maps are differenced to estimate the map noise. The 
95 GHz and 150 GHz full-depth maps of a given field 
are then combined using a spatial-spectral matched fil- 
ter (e.g., Melin et al.|2006 l that optimizes signal-to-noise 
on cluster-shaped objects with an SZ spectral signature. 
Cluster candidates are identified in the resulting filtered 
map using a simple peak detection algorithm, and each 
candidate is assigned a signal-to-noise value based on the 
peak amplitude divided by the RMS of the filtered map in 
the neighborhood of the peak. Twelve different matched 
filters are used, each assuming a different scale radius 
for the cluster, and the maximum signal-to-noise for a 
given candidate across all filter scales is referred to as ^, 
which we use as our primary SZ observable. In 2008, the 
95 GHz detectors in the SPT receiver had significantly 
lower sensitivity than the 150 GHz array, and the cluster 
candidates from those observations are identified using 
150 GHz data only; the candidates from 2009 observa- 
tions were identified using data from both bands. The 
data from the two observing seasons yielded a total of 
224 cluster candidates with ^>4.5 — the sample discussed 
here. 

2.2. Optical/NIR Imaging 

The cluster candidates detected using the method de- 
scribed above are followed up by optical and, in many 
cases, NIR instruments. In this section, we describe 
the overall optical/NIR follow-up strategy, the differ- 
ent imaging and spectroscopic observations and facilities 
used, and the data reduction methods used to process 
the raw images to catalogs. 

2.2.1. Imaging Observations 

The optical/NIR follow-up strategy has evolved since 
the first SPT-SZ candidates were identified. Originally 
we imaged regions of the sky with uniformly deep, multi- 
band observations in griz optical bands to confirm SZ de- 
tections and estimate redshifts as in Staniszewski et al.l 
(j2009). For the first SPT cluster candidates, we used 
imaging from th e Blanco Cosmology Survey (BCS; De 
Isai et al.l 120121) to follow up candidates in parts ot 



the 2008 fields. The BCS is a 60-night, -80 deg^ 
NOAO survey program carried out in 2005-2008 using 
the Blanco/MOSAIC2 griz filters. The BCS survey was 
completed to the required depths for 5a detection at 
OAL* within 2'.'3 apertures up to z— 1. The goal of this 
survey was to provide optical imaging over a limited area 
of the SPT survey to enable rapid optical follow-up of the 
initial SPT survey fields. 

For clusters outside the BCS region we initially ob- 
tained deep griz imaging on a cluster-by-cluster basis. 
But as the SPT survey proceeded and the cluster candi- 
date list grew, it became clear that this strategy was too 



costly, given the limited access to follow-up time. More- 
over, eventually the full SPT region will be imaged to 
uniform lOcr depths of mag~24 in griz by the DES. We 
therefore switched to an adaptive strategy of follow-up 
in which we observed each SPT cluster candidate to the 
depth required to find an optical counterpart and deter- 
mine its redshift. 

For each SPT cluster candidate, we perform an ini- 
tial pre-screening of candidates using the Digitized Sky 
Survey (DSS)"^^. We examine DSS images for each clus- 
ter candidate to determine whether it is "low-z" or 
"higli-z" , where the redshift boundary lies roughly at 
z = 0.5. We use the DSS designation to prioritize the 
target list for the appropriate telescope, instrument and 
filters with which we observe each candidate. Specif- 
ically, candidates that are clearly identified in DSS im- 
ages are likely to be low-z clusters and are designated for 
follow-up observations on small-aperture (lm-2m) tele- 
scopes, while high-z candidates are designated for large- 
aperture (4m-6.5m) telescopes. The various ground- and 
space-based facilities used to collect optical/NIR imaging 
data on SPT clusters are summarized in Table [H Each 
telescope/instrument combination is assigned a numeric 
alias that is used to identify the source of the redshift 
data for each cluster in Table [3l 

For the > 4m-class observations, we use an adaptive 
filter and exposure time strategy so that we can efficiently 
bracket the cluster member galaxy's 4000A break to the 
depth required for redshift estimation. In this approach 
we start with a first imaging pass, where each candidate 
is observed in the g, r, and z bands to achieve a depth 
corresponding to a 5a detection of a 0.4L* galaxy at 
z — 0.8 in r and z bands. Observations are also taken in 
a single blue filter for pho tomet ric calibration using the 
stellar locus (discussed in { 2.2.2 1. For candidates with no 
obvious optical counterpart after first-pass observations, 
a second-pass is executed to get to z = 0.9 in i and z 
bands. 

If the candidate is still not confirmed after the second- 
pass in i and z bands, and is not covered by the 
Spitzer/IRAC pointed observations described below, we 
attempt to obtain ground-based NIR imaging for that 
candidate using the NEWFIRM camera on the CTIO 
Blanco telescope. The data presented here are imaged 
with NEWFIRM during three observing runs in 2010 
and 2011, yielding Kg data for a total of 31 candidates. 
Typical observations in Ks consist of 16 point dither pat- 
terns, with 60 second exposures divided among 6 coadds 
at each dither position. Median seeing during the 2010 
runs was 1'.'05; during the 2011 run observing conditions 
were highly variable and the seeing ranged from 1'.'05 to 
2'.'6. 

We note that most of the galaxy cluster candidates 
in this work wit h significance £ > 4.8 were imaged 
with Spitzer (We rner et a l.l 2004). More specifically, 
Spitzer/IRAC ( jFazio et ai. 2004 ) imaging has been ob- 
tained for 99 SZ cluster candidates in this work. The 
on-target Spitzer observations consist of 8 x 100 s and 
6x30 s dithered exposures at 3.6/im and 4.5/im, respec- 
tively. The deep 3.6/Ltm observations should produce 5a 
detections of passively-evolving O.IL* cluster galaxies at 
z = 1.5. 

http://archive.stsci.edu/dss/ 
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For some of the NIR analysis, we augment the data 
from our Spitzer and NEWFIRM observations with the 
recently releas ed all-sky Wide-fie ld Infrared Survey Ex- 
plorer ( WISE; [Wright et al.|2010[ ) data. Finally, we note 
that a few of the clusters were observed with Magel- 
lan/Megacam t o obtain weak gra vitational lensing mass 
measurements (High et al. 2012). These data are natu- 
rally much deeper than our initial followup imaging. 

2.2.2. Data Processing 

We use two independent optical data processing sys- 
tems. One system, which we refer to as the PHOTPIPE 
pipeline, is used to process all optical data except Mag- 
ellan/Megacam data, and the other, which is a devel- 
opment version of the Dark Energy Survey data man- 
agement (DESDM) system, is used only to process the 
Blanco/Mosaic2 data. PHOTPIPE was used t o process 
optical data for previous SPT cluster catalogs dVander- 



linde et al |20Tol|High et aL|2010[[Wimamson et al.|201l[ ); 
the DESDM system has been used as a cross-check in 
th ese works and was the prim ary reduction pipeline used 



Staniszewski et al. 
he basic stages of 



(|2009). 

HeTHOTPIPE pipeline, initially 



developed for the SuperMACHO an d ESSENCE project s 



and described in Rest et al. (2005), Garg et al. 



(2007), 



and Miknaitis et al.| ( |2007[ ), include tlat-helding, astrom- 
etry, coadding, and source extraction. F urther details 
are g iven in HIO. In the DESDM system (|Ngeow et al.| 
2006| pohr ct al.'2008; Dcsai ct al. 2012), the data from 
each night tirst undergo detrending corrections, which 
include cross-talk correction, overscan correction, trim- 
ming, bias subtraction, flat fielding and illumination cor- 
rection. Single epoch images are not remapped to avoid 
correlating noise, and so we also perform a pixel-scale 
correction that brings all sources on an image to a com- 
mon photometric zeropoint. For i and z bands we also 
carry out a fringe correction. Astrometric calibration is 
done by using the AstrOmatic code SCAMP ( Bertin|2006 1 
and the USNO-B catalog. Color terms to transform to 
the SDSS system rely on photometric solutions derived 
from observations of SDSS equatorial fields during pho- 



tometric ni 



ghts (jDe 



et al. 



In both pipeline s 



2012). 



coaddition 
the 



is done using 
DESDM system 



SWarp ( [Bertin et al.||2002[ ). In 

the single epoch images contributing to the coadd are 
brought to a common zeropoint using stellar sources 
common to pairs of images. The final photometric 
calibration of the coadd images is carried out using the 
stellar locus as a constraint on t he zeropoint offset s 



2009), 



between neighboring bands (e.g., [High et al. 
while the absolute photometric calibration comes from 
2MASS (Skrutskie ct al. 2006) ). For griz photometry the 
calibration is done with reference to the median SDSS 



stellar locus ( jCovey et al. 2007), but for the Swope 
data using Johnson filters, the calibration relies on a 
stellar locus derived from a se quence of models of stellar 
atmospheres from PHOENIX ( Brott fc Hauschildt|20 05 1 
with empirically measured CCd, filter, and atmosphere 
responses. Catalo ging is done using SExtractor (Bertin, 



fc Arnouts 
we separate. 



1996), 



and within the DESDM catalogs 

y calibrate mag-auto and a 3" aperture 

magnitude using the stellar locus. The small aperture 

photometry is used to obtain higher signal to noise 
galaxy color measurements. 



Quality checks of the photometry are carried out on a 
cluster by cluster basis using the scatter of stars about 
the expected stellar locus and the distribution of off- 
sets in the single-epoch photometry as a function of cal- 
ibrated magnitude (so-called photometric repeatability 
tests). Poor quality data or failed calibrations are easily 
identified as those coadds with high stellar locus scatter 
and or high scatter in the photometric repeatability tests 
(see lDesa et al.||2012[ ) . 

NEWI'IRM imag ing data are reduce d using the FAT- 
BOY pipeline (Eikenberry et al.||2006|), originally devel- 
oped for the FLAMINCOS-2 instrument, and modified 
to work with NEWFIRM data in support of the Infrared 



Bootes Imaging Survey (Gonzalez et al. 2010). Indi 



vidual processed frames are combined using SCAMP and 
SWarp, and photometry is calibrated to 2MASS. 

Spitzer/IRAC data a re reduced following the method 
of Ashby et al. (20091. Briefiy, we correct for column 
pulldown, mosaic the individual exposures, resample to 
(y.'86 pixels (half the solid angle of the native IRAC pix- 
els), and reject cosmic rays. Magnitudes are measured in 
4"~diameter apertures and corrected for the « 40% loss 
due to the broad PSF. The Spitzer photometry is crucial 
to the measurement of photometric redshifts for clusters 



at z > 0.8, as described in ^3.1 



The acquisition and processing for the initi al weak lens- 
ing M egacam data is described in detail in 'High et"aL 
(2012 1. These data are reduced separately from the other 



imaging data using the Smithsonian Astrophysical Ob- 
servatory (SAO) Megacam reduction pipeline. Standard 
raw CCD image processing, cosmic-ray removal, and flat- 
fielding are performed, as well as an additional illumina- 
tion correction to account for a low-order scattered light 
pattern. The final images are coadded onto a single pixel 
grid with a pixel scale of 0'.'16 using SWarp. Sources are 
detected in the coadded data in dual-image mode using 
SExtractor, where the r band image serves as the de- 
tection image. The photometry is calibrated by fitting 
colors to the stellar locus, and color-term corrections are 
accounted for in this step. 



2.3. Spectroscopic Observations 

We have targeted many of the SPT clusters with long- 
slit and multi-object spectroscopy, and some of the spec- 
troscopic redshifts have appeared in previous SPT publi- 
cations. We have used a variety of instruments: GMOS- 
S^"^ on Gemini South, F0RS2 ( [Appenzeller et al.||l998 l 
on VLT, LDSS3 on Magellan- Clay, and the IMACS cam- 
era on Magellan Baade (in long-slit mode and with the 
GISMO^*^ complement). 

A detailed description of the configurations, observ- 
ing runs and reductions will be presented elsewhere (J. 
Ruel et al., in prep). In all cases, for a given cluster 
we target bright galaxies measured to have a high prob- 
ability of lying on the cluster's red sequence and ob- 
served these galaxies with a combination of filter and 
disperser that yielded a low-resolution spectrum around 
their Ca II H&K lines and 4000A break. CCD reduc- 
tions are made using standard packages, including COS- 



http:/ /www. gemini.edu/node/10625 
http:/ /www.lco.cl/telescopes-information/ 
magellan / instruments / imacs /gismo /gismoquickmanual.pdf 
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TABLE 1 
Optical and infrared imagers 



Ref.'' 


Site 


Telescope 


Aperture 
(m) 


Camera 


Filters'' 


Field 


Pixel scale 
(") 


1 


Cerro Tololo 


Blanco 


4 


MOSAIC-II 


griz 


36' X 36' 


0.27 


2 


Las Campanas 


Magellan/Baade 


6.5 


IMACS f/2 


griz 


27^4 X 27.'4 


0.200 


3 


Las Campanas 


Magellan/Clay 


6.5 


LDSS3 


griz 


8f3 diam. circle 


0.189 


4c 


Las Campanas 


Magellan/Clay 


6.5 


Megacam 


gri 


25' X 25' 


0.16 


5 


Las Campanas 


Swope 


1 


SITeS 


BVRI 


14f8 X 22.'8 


0.435 


6 


Cerro Tololo 


Blanco 


4 


NEWFIRM 


Ks 


28' X 28' 


0.4 


7 




Spitzer Space Telescope 


0.85 


IRAC 


3.6/im, 4.5/im 


5f2 X 5f2 


1.2 


8 




Wide-field Infrared Survey Explorer 


0.40 




3.4/im, 4.6^m 


47' X 47' 


6 



Note. — Optical and infrared cameras used in SPT follow-up observations. 

Shorthand alias used in Tabic [3] 

Not all filters were used on every cluster. 

Megacam data were acquired for a large follow-up wcak-lensing program. 



MOS ( |Kelson| 12003) for IMACS data and IRAF^^ for 
CMOS and~F0RS2. Redshift measurements are made 



by cross-cor relation with the RVSAO package (Kurtz & 
Mink 1998) and a proprietary template fitting method 
that uses SDSS DR2 templates. Results are then visu- 
ally confirmed using strong spectral features. 

In Table [2] the source for every spectroscopic redshift 
is listed, along with the number of cluster members used 
in deriving the redshift. For clusters for which we report 
our own spectroscopic measurements, we list an instru- 
ment name and observation date; we give a literature 
reference for those for which we report a value from the 
literature. In Tabic [3j we report spectroscopic redshifts 
for 57 clusters, of which 36 had no previous spectroscopic 
redshift in the literature. Unless otherwise noted, the 
reported cluster redshift is the robust biweight average 
of the redshifts of all spectroscopically confirmed cluster 
members, and the cluster redshift uncertainty is found 
from bootstrap resampling. 



3. METHODOLOGY 

In this section, we describe the analysis methods used 
to: 1) extract cluster redshift estimates and place red- 
shift limits; 2) empirically verify the estimates of catalog 
purity; and 3) measure rBCG positions. 



the magnitudes of similar galaxies. In the optical a nal- 
ysis for our two previous cluster catalog releases ( |High| 
al.|2010[|Williamson et al.|2011[ ), we relied on MetEod 
tor the results and Method 3 as a cross-check. In this 
work, we improve the precision of the measured redshifts 
by applying multiple redshift estimation algorithms and 
combining the results. Through cross-checks during the 
analysis, we find that these methods have different fail- 
ure modes and that comparing the results provides a way 
of identifying systems that require additional attention 
(including systems where the cluster's central region is 
contaminated by foreground stars or the cluster resides 
in a crowded field). 

All thr ee methods use th e singl e-stell ar-population 
models of [Bruzual fc Charlotj (BC03; |2003[ ). These mod- 
els allow us to transform the location of the red-sequence 
overdensity in color space to a redshift estimate. A model 
for the red galaxy population as a function of redshift is 
built assuming a single burst of star formation at = 3 
followed by passive evolution thereafter. Models are se- 
lected over a range of metallicities and then calibrated 
to reproduce the color and tilt of the re d sequence in 
the Coma cluster (Eisenhardt et al. 2007) at z = 0.023. 



The calibra t ion p rocedure is described in more detail in 



Song et al. (2012). These models are used in determin- 



ing exposure times and appropriate filter combinations 
for imaging observations, and in the calculation of red- 
shifts and redshift limits from those observations. 



3.1. Photometric Redshifts 

Using the procedure described in §2.2[ we obtain 
ground-based imaging data and galaxy catalogs that in 
most cases allow us to identify an obvious overdensity of 
red-sequence galaxies within approximately an arcminute 
of the SPT candidate position. For these optically con- 
firmed cluster candidates, we proceed to estimate a pho- 
tometric redshift. 

In this work, we employ three methods (which we refer 
to Method 1, 2, and 3 in the following sections) to esti- 
mate cluster redshifts from optical imaging data. Two 
methods (Method 1 and 2) use the color of the galaxies in 
the cluster red sequence, and the third (Method 3) uses 
the average of red-sequence galaxy photometric redshifts 
estimated with a neural-network algorithm, trained with 

http://iraf.noao.edu 



3.1.1. Photometric Redshift Measurement Methods 

In Method 1, a cluster is confirmed by identifying an 
excess of galaxies with colors consistent with those de- 
rived from BC03 (simultaneously for all observed filters), 
after subtracting the background surface density. The 
background-subtracted galaxy number is extracted from 
an aperture within a radius of (3.5,2.5,1.5)' from the SPT 
candidate position and uses galaxies with photometric 
color uncertainties < (0.25,0.35,0.45) and apparent mag- 
nitudes brighter than m*-|-(3,2,l) (or the magnitude limit 
of the data) in the red sequence, for z<0.2, 0.2<z< 0.6, 
and z>0.6 respectively. The background measurement 
is obtained by applying the same criteria outside of the 
cluster search aperture. The redshift is estimated from 
the most significant peak in this red-sequence galaxy ex- 
cess. Improvements over the implementation in HIO in- 
clude using additional colors (r-z and g-i, plus NIR col- 
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TABLE 2 
Spectroscopic Follow-Up 





Cluster 




Inst 


Obs 


# 


Refs 




Cluster 




Inst 


Obs 


# 


Refs 


SPT 


CLJOOOO 


5748 




oCp 


om n 

zUlU 


26 




SPT 


CLJ2058 


5608 


T\ T/^O o 

GMOb-b 


Sep 2011 


9 




SPT 


CLJ0205 


5829 


T TV T A O 

IMACb 


oep 




9 


a 


SPT 


CLJ2100 


4548 


T~l/^ T^ O O 

FORS2 


Aug 2011 


19 




SPT 


CLJ0205 


6432 


GMOb-b 


oep 




15 




SPT 


CLJ2104 


5224 


lORbi 


Jun 2011 


23 




SPT 




5819 


GMOb-s 


oGp 




10 




SPT 


T T o "1 
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5844 


T71/~\T^ O O 

lORbi 


Dec 2010 


15 




SPT 


CLJ0234 


5831 


TO TV T 


wct 




22 


b 








TO TV T 

GIbMO 
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3 


i 


SPT 


CLJ0240 


5946 


TO TV 




901 n 


25 




SPT 


CLJ2118 


5055 


lORbi 


TV /T or\ 1 1 

May 2011 


25 




SPT 


CLJ0254 


5857 


TOTV n r\ 


Oof 

wet 




35 


b 


SPT 


CLJ21i4 


6124 


TO TV n r\ 

GIbMO 


Sep 2009 


24 




SPT 


CLJ0257 


5732 


TO TV If r\ 


wct 




22 




SPT 


CLJ2130 


6458 


TOTV n r\ 

GIbMO 


bep 2009 


47 




SPT 


CLJ0317 


5935 


TO TV T 

GIbMO 


Oof 




17 




SPT 


CLJ2135 


5726 


TO TV T /~\ 

GIbMO 


Sep 2010 


33 




SPT 


CLJ03i8 


5541 










c 


SPT 


T T o "1 n 

CLJ2136 


4704 


T\ T/~\0 o 

GMOb-b 


Sep 2011 


24 




SPT 


CLJ0431 


6126 










c 


SPT 


CLJ2136 


6307 


TO TV If r\ 

GIbMO 


Aug 2010 


10 




SPT 


CLJ0433 


5630 


TO TV T 

GIbMO 


J an 


om 1 


22 




SPT 


CLJ2138 


6007 


TO TV T /~\ 

GIbMO 


Sep 2010 


34 




SPT 


CLJ0509 


5342 


<^ T\ fT/~\0 O 

GMUb-b 


uec 


zuuy 


18 


d, e 


SPT 


<^T TO -1 AC 

CLJ2145 


5644 


TO TV T/^ 

GIbMO 


bep 2009 


37 




SPT 


CLJ0511 


5154 


<^ T\ T/~\0 O 

GMOb-b 


oep 


om 1 


15 




SPT 


CLJ2146 


4633 


TTV If \ /^O 

IMACb 


Sep 2011 


17 




SPT 


CLJ0516 


5430 


TO TV T 

GIbMO 


oep 




48 


f 


SPT 


T T O -1 ^ O 

CLJ2146 


4846 


TV T O O 

GMOS-S 


Sep 2011 


26 




SPT 


CLJOoil 


5104 










e 


SPT 


T T O "1 /I O 

CLJ2148 


6116 
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GIbMO 


bep 2009 


30 




SPT 


CLJOoio 


5300 


T\ T/~\o o 

GMOb-b 


J an 


Zulu 


20 


d, e 


SPT 


CLJ2155 


6048 


*^ T\ T/^O O 

GMOb-b 


Sep 2011 


25 




SPT 


CLJ0533 


5005 
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4 


d 


SPT 
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CLJ2201 


5956 








c 


SPT 


CLJ0534 


5937 


LDSS3 
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2008 


3 




SPT 


CLJ2300 


5331 
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GIbMO 


Oct 2010 


24 




SPT 


CLJ0546 


5345 
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21 


cr 
O 


SPT 


CLJ2301 


5546 


GISMO 


Aug 2010 


11 










GMOS-S 
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2 


e 
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CLJ0551 
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34 
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50 


d 


SPT 


CLJ0559 
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37 


d, e 


SPT 


CLJ2332 


5358 


GISMO 
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24 




SPT 


CLJ2011 


5725 










f 


SPT 


CLJ2337 


5942 


GMOS-S 


Aug 2010 


19 


d 


SPT 


CLJ2012 


5649 










c 


SPT 


CLJ2341 


5119 


GMOS-S 


Aug 2010 


15 


d 


SPT 


CLJ2022 


6323 


GISMO 


Oct 


2010 


37 




SPT 


CLJ2342 


5411 


GMOS-S 


Sep 2010 


11 




SPT 


CLJ2023 


5535 










f 


SPT 


CLJ2351 


5452 








h 


SPT 


CLJ2032 


5627 


GISMO 
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2010 


31 




SPT 


CLJ2355 


5056 


GISMO 


Sep 2010 


37 




SPT 


CLJ2040 


5725 


GISMO 
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2010 


5 




SPT 


CLJ2359 


5009 


GISMO 


Aug 2010 


21 




SPT 


CLJ2043 


5035 


FORS2 


Aug 


2011 


21 










GMOS-S 


Dec 2009 


9 




SPT 


CLJ2056 


5459 


GISMO 


Aug 


2010 


12 



















Note. — Instruments [Inst] : GMOS-S on Gemini South 8m, IMACS on Magellan Baade 6.5m, GISMO complement to 
IMACS on Magellan Baade 6.5m, LDSS3 on Magellan Clay 6.5, FORS2 on VLT Antu 8m; Observing date s [Obs]: dates 
each data taken; Number of galax ie s [#]: Number of galaxie s used in deriv ing re dshifts; Referen ces [R efs ]: °|Stalder et al. 

" Williamson et al.H201l|l, jStruble fc Rood] jl999l , lHigh et~ l|2010t , ' [Sifon et ab] l |2012[ , j |Bohringer et al, 
Br odwin et al.| l |2010[ l, " |Buckley-Geer et a'L7 l |201l( , 1 Foley et ab] " 201l[ l 



2012 



2004 



ors) in the red-sequence fitting, using the deeper photom- 
etry available from coadded images, and sampling the 
entire CCD mosaic rather than a single CCD for better 
background estimation. 

Method 2 is similar in that it searches for an over- 
density of red-sequence galaxies. This method, used to 
estimate the redshifts for a sample of 46 X-ray-selected 
clusters ( Suhada et al. 2012 ) is d escribed and tested in 
more detail in Song et al. (2012 1. It includes a mea- 
sure of the background surtace density based on the 
entire imaged sky area surrounding each cluster candi- 
date and subtracts the background from the red galaxy 
counts in an aperture of 0.8 Mpc. Only galaxies with 
luminosity > 0.4L* and magnitude uncertainty < 0.25 
are used, and the aperture and luminosity are recalcu- 
lat ed for each potenti al redshift. Originally as described 
in ISong et al. (2012), we search for an overdensity of 



red-sequence galaxies using two or three available color- 
redshift combinations simultaneously for every cluster; 
essentially, we scan outwards in redshift using the follow- 
ing color combinations: g-r and g-i for z <0.35, g-i and 
r-i and r-z for 0.3<z<0.75, r-z and i-z for z>0.75. The 
cluster photometric redshift is extracted from the peak of 
the galaxy overdensity in redshift space. The redshift is 
then refined by fitting the red-sequence overdensity dis- 
tribution in redshift space with a Gaussian function. The 
ve rsion used here (wh ich is the same as the method used 
in Suhada et al.|20l"2 | has one more refinement, in which 



the colors of the galaxies that lie in the peak redshift bin 
identified by the overdensity method are converted into 
individual galaxy photometric redshifts. In this conver- 
sion we assume that the galaxies are red-sequence cluster 
member galaxies, and the photometric redshift uncer- 
tainty reflects the individual photometric color errors. A 
final cluster redshift is calculated as an inverse- variance- 
weighted mean of these galaxy photometric redshifts. 

Method 3 shares the same principle as the other two 
in that it involves searching for a density peak in the 
galaxy distribution near the position of the SPT candi- 
date. We first select individual red cluster members us- 
ing location relative to the SPT candidate position and 
galaxy color as the criteria for cluster membership. For 
the redshifts presented here, this is done visually using 
pseudo-color coadded images for each cluster, although 
in principle this could be automated. Galaxy selection 
is not confined by a specific radial distance from SZ cen- 
ters as in the other two methods, nor by photometric 
uncertainties. Selected galaxie s are then fe d into Artifi- 
cial Neural Networks (ANNz; 'CoUi ster fc~Lahav||2004| ), 
which is trained using the same BC03 models used in 
the other methods. ANNz returns redshift estimates for 
individual galaxies, and a peak in galaxy redshift distri- 
bution is adopted as the initial cluster redshift. Then, 
as in Method 2, individual galaxy photometric redshifts 
are averaged using inverse- variance weighting to produce 
the cluster photometric redshift. With this initial esti- 
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mate of the redshift, we then perform an outher rejection 
using iterative la chpping, where the la corresponds to 
the root-mean-square (RMS) variation of the measured 
galaxy photometric redshift distribution. Once the re- 
jection is carried out, we refine the cluster photometric 
redshift estimate using the weighted mean of the non- 
rejected sample of cluster galaxies. No outlier rejection 
is undertaken if there are fewer than 20 selected galaxies 
in the original sample. 

Method 3 is a good cross-check, as well as a stand- 
alone redshift estimator, because we can visually confirm 
which galaxies contribute to the redshift determination. 
Although this method requires photometry in more than 
just two bands, it appears to be less susceptible to the 
problems in two-band methods that are associated with 
pileup of red sequence galaxies at redshifts where the 
4000A break is transitioning out of a band. 
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Fig. 1. — Top: Photometric redshift 2phot versus spectroscopic 
redshift 2apcc for each redshift estimation method for 47 spectro- 
scopicaUy confirmed clusters at z < 0.9 where we use only griz 
photometry. Bottom: the distribution of the photometric redshift 



residuals Az 



2phot 



as a function of 



^spcc ■ 



Inset: the nor- 



malized residual distributions, which all have RMS{Az/az^j^^^) ~ 
1. The RMS scatter of Az/{1 + z) is 0.028, 0.023 and 0.024 for 
Methods 1, 2 and 3, respectively. 



Next we characterize redshift estimates from each 
method using spectroscopically confirmed clusters. We 
use 47 clusters with spectroscopic redshifts (zspcc) where 
only griz data are used for photometric redshift estima- 
tion. In this process, photometric redshift (zphot) biases 
(namely, smooth trends of photometric redshift offset as 
a function of redshift) are measured and corrected in 
Method 1 and 2, while no significant bias correction is 
necessary for Method 3. Once biases are removed, we 



examine the photometric-to-spectroscopic redshift off- 
sets to characterize the performance of each method. 
We find the RMS in the quantity A2:/(1-|-Zspcc), where 
Az = Zphot — ^^spec to have values of 0.028, 0.023, and 
0.024 in Methods 1, 2, and 3, respectively (see Figure[T]). 

Our goal is not only to estimate accurate and precise 
cluster redshifts, but also to accurately characterize the 
uncertainty in these estimates. To this end, we use the 
spectroscopic subsample of clusters to estimate a system- 
atic floor (Tgys in addition to the statistical component. 
We do this by requiring that the reduced describing 
the normalized photometric redshift deviations from the 
true redshifts Xrcd ~ 

X;(Az/cr2p,^^j2/7Vdof have a value 
Xicd ~ 1 fo'^ each method, where az^t.^^ is the uncertainty 
in measured Zphot and iVdof is the number of degrees of 
freedom. We adopt uncertainties a'^^^^^ = a^^^^. + a'^^^ 

and adjust asys to obtain the correct Xred- ^^^^ tuning 
process we also include redshift estimates of the same 
cluster from multiple instruments when that cluster has 
been observed multiple times. This allows us to test the 
performance of our uncertainties over a broader range of 
observing modes and depths than is possible if we just 
use the best available data for each cluster. 

For Method 1 we separately measure the systematic 
floor cTgys for each different photometric band set. For 
the grizKs instruments (Megacam, IMACS, LDSS3, 
M0SAIC2, NEWFIRM), we estimate a^ys = 0.039 
; for the BVRI instrument (Swope), a^ys ~ 0.033; 
and for Spitzer-oivXy, asys = 0.070. In Method 2, we 
find (jgys = 0.030 for the griz instruments (Megacam, 
IMACS, LDSS3, M0SAIC2). For Method 3 we estimate 
fgys = 0.028 for the griz instruments. 

Once this individual estimation and calibration is done, 
we conduct an additional test on the redshift estimation 
methods, again using the spectroscopic subsample. The 
purpose of this test is to see how the quality of pho- 
tometry (i.e., follow-up depth) affects the estimations. 
We divide the spectroscopic sample into two groups: in 
one group, the photometric data is kept at full depth, 
while the photometric data in the other group is manu- 
ally degraded to resemble the data from the shallowest 
observations in the total follow-up sample. To create the 
'shallow' catalogs, we add white noise to the full-depth 
coadds and then extract and calibrate catalogs from these 
artificially noisier images. Results of this test show that 
the accuracy of the photometric redshift estimation is 
affected by the poorer photometry, but that this trend 
is well captured by the statistical uncertainties in each 
estimation method. 

3.1.2. Combining Photo-z Estimates to Obtain Zcomb 

Once redshifts and redshift uncertainties are estimated 
with each method independently, we compare the differ- 
ent redshift estimates of the same cluster. Note that this 
comparison is not possible for Swope or Spitzer-only red- 
shifts, which are measured only with Method 1. Outliers 
at > 3cr (>6%) in H-Zphot are identified for additional 
inspection. In some cases, there is an easily identifiable 
and correctable issue with one of the methods, such as 
misidentification of cluster members. If, however, it is 
not possible to identify an obvious problem, the outliers 
are excluded from the combining procedure. This out- 
lier rejection, which occurs only in two cases, causes less 
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than 0.05 change in the combined Zphotin both cases. 

We combine the individual estimates into a final best 
redshift estimate, Zcomb, using inverse- variance weighting 
and accounting for the covariance between the methods, 
which we expect to be non-zero given the similarities in 
the methods and the common data used. Correlation co- 
efBcients for the photometric redshift errors among the 
different methods are measured using the spectroscopic 
sample. The measured correlation coefficient, r^, be- 
tween each pair of methods is 0.11 (Method 1 & 2), 0.40 
(Method 2 & 3) and 0.19 (Method 1 & 3). 

With the correlation coefficients we construct the op- 
timal combination of the individual estimates as: 



1 



^comb 



(1) 



where Wij — C^j ^ , and the covariance matrix Cij is com- 
prised of the square of the individual uncertainties along 
the diagonal elements (af) and the product of the mea- 
sured correlation coefficient and the two individual un- 
certainty components (i.e., TijaiUj) on the off-diagonal 
elements. The associated uncertainty is 



1 



(2) 



Because of the positive correlations between the three 
methods' errors, the errors on Zcomb are larger than would 
be the case for combining three independent estimates; 
however, we do see an improvement in the performance 
of the combined redshifts relative to the individual es- 
timates that is consistent with the expectation given 
the correlations. The performance of this combined red- 
shift method is presented in Figure |2] the residual dis- 
tribution is roughly Gaussian, and the associated uncer- 
tainties provide a good description of the scatter of the 
redshift estimates about the spectroscopic redshifts (the 
RMS variation of Az/cr^pj^^^ is 1.04). The benefit from 
combining different measurements is evidenced from the 
tighter distribution in the redshift versus Zgpec plot; the 
RMS scatter of Az/(1 -f Zspoc) is 0.017, corresponding 
to a ~40% improvement in the accuracy relative to the 
accuracy of a single method. 

3.1.3. Spitzer Photometnc Redshifts 

For clusters where we do not have deep enough optical 
data to estimate a redshift but that do have Spitzer cov- 
erage, we use Method 1 to measure the redshifts using 
Spitzer-only colors in the same manner as we do with op- 
tical data. Overdensities of red-sequence galaxies in clus- 
ters have been identified using Spitzer-only color selection 
at high redshift, where the IRA C bands are probing the 
peak of the stellar emission (Stern et al.||2005" Papovich 
20081. In our sample, the comparison of Spitzer-only 
redshifts with spectroscopically derived redshifts shows 
good performance, indicating that the assumption of a 
wcU-dcvcl oped red- sequcncc app ears to hold out to z > 1 



(e.g.. Bower et al. 



1992, Eisenhardt et al. 2008 



Muzzin 



al. 2009). INote that the possibility ot the cluster bemg 
al lower redshift is already ruled out from the available 
optical data for these candidates. 

Figure [3] shows the performance of the Spitzer-only 
redshifts in eight clusters where spectroscopic redshifts 




0.2 0.4 0.6 

Spectroscopic redshift 

Fig. 2. — Top: Weighted mean photometric redshift Zcomb versus 
spectroscopic redshift using the same sub-sample as in Figure ^ 
Bottom: the distribution of the redshift errors. The RMS scatter m 
A2/(l-|- Zspoc)=0.017. Inset: histogram of the normahzed redshift 
error distribution, which is roughly Gaussian with RMS ~ 1. 



are available. Although the accuracy in Zphot is lower 
(Az/(1 -\- z) 0.049) than those derived from optical- 
only or optical-IRAC colors, the performance is reliable. 
We fiag these cases in the final table to make note of this 
difference in method. The larger uncertainties of Spitzer- 
only derived redshifts are possibly due to the broad width 
of IRAC filters and the fact that AGN emission or vig- 
orous star formation can shift the location of the l.dfini 
bump. 

3.1.4. Redstiift Limits 

In most cases there is an obvious, rich overdensity of 
red cluster galaxies in our follow-up imaging, from which 
it is straightforward to confirm the galaxy cluster giving 
rise to the SZ signal and to estimate the cluster redshift. 
For unconfirmed candidates, it is not possible to say with 
absolute certainty that no optical/NIR counterpart ex- 
ists; with real, finite-depth optical and NIR data, the 
possibility always exists that the cluster is distant enough 
that no counterpart would have been detected at the 
achieved optical/NIR depth. Assigning a relative prob- 
ability to these two interpretations of an optical/NIR 
non-detection (i.e., a false detection in the SZ data or a 
higher-redshift cluster than the optical/NIR observations 
could detect) is especially important for interpreting the 
SZ cluster sample cosmologically. To this end, we calcu- 
late a lower redshift limit for every SZ-selected candidate 
for which no counterpart has been found. Because the 
optical/NIR follow-up data is not homogeneous, we do 
this separately for each unconfirmed candidate. 

To estimate the depth of the optical/NIR coadded im- 
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Fig. 3. — Photometric redshift vs. Zspcc ior clusters where only 
the Spitzer IRAC 3.6/im-4.5/.tm colors are used. In all cases where 
we present Spitzer photo-z's, we have optical data to rule out the 
presence of a low-redshift cluster. 



ag es, we utilize a Mon te-Carlo based technique described 



Ashby et al. ( 2009 ) . In brief, we perform photometry 
of the sky in various apertures at 1500 random positions 
in each image. To measure the sky noise, we then fit a 
Gaussian function to the resulting flux distribution (ex- 
cluding the bright tail which is biased by real sources 
in the image). Taking the measured sky noise from 3"- 
diameter apertures, we add a PSF-dependent aperture 
correction. A redshift limit is derived for each filter by 
matching a 0.4L* red-sequence galaxy from the BC03 
model to the measured lOcr magnitude limit. We use the 
redshift limit from the second deepest filter with regards 
to 0.4L* red-sequence objects, as we require a minimum 
of two filters to measure a redshift. These redshift limits 
are compared to limits derived by comparing observed 
number counts of galaxies as a function of mag nitude to 
distribution s derived from much deeper data (Zenteno 
et al.|[20TT ). We find the two independent redshitt iimit 
estimations are in good agreement. 

For cluster candidates with 5piteer/IRAC observa- 
tions, the redshift estimation is not limited by the depth 
of the optical data, and we use the IRAC data to cal- 
culate a lower redshift limit for these candidates. The 
IRAC data are highly uniform, with depth sufficient to 
extract robust photometry down to Q.IL* out to a red- 
shift oi z = 1.5. In principle, ~ 0.5L* photometry should 
be sufficient for redshift estimation; however, we adopt 
z = 1.5 as a conservative lower redshift limit for any 
unconfirmed candidates with IRAC data. 

3.2. NIR Overdensity Estimates for Unconfirmed 
Candidates 

For cluster candidates for which we are unable to esti- 
mate a redshift, we can in principle go beyond a simple 
binary statement of "confirmed/unconfirmed" using NIR 
data. Even if there is not a sufficient number of galaxies 
in the NIR data to estimate a red sequence, there is in- 
formation in the simple overdensity of objects (identified 
in a single NIR band) within a certain radius of an SPT 
candidate, and we can use this information to estimate 
the probability of that candidate being a real, massive 



cluster. We can then use this estimate to sharpen our 
estimate of the purity of the SPT-selected cluster sam- 
ple. We calculate the single-band NIR overdensity for 
all unconfirmed candidates using WISE data, and we 
compare that value to the same statistic estimated on 
blank-field data. We perform the same procedure using 
Spitzer/IRAC and NEWFIRM data for unconfirmed can- 
didates that were targeted with those instruments. For 
comparison, we also calculate the same set of statistics 
for each confirmed cluster above z = 0.7. 

We estimate the galaxy overdensity within a 1 arcmin^ 
aperture. To increase the signal-to-noise of the estima- 
tor, we assume an angular profile shape for the cluster 
galaxy distribution and fit the observed distribution to 
this shape. The assumed galaxy density profile is a pro- 
jected /3 model with /3 = 1 (the same profile assumed for 
the SZ signal in the matched-filter cluster detection al- 
gorithm in R12). We have tried using a projected NFW 
profile as well, and the results do not change in any sig- 
nificant way (due to the relatively low signal-to-noise in 
the NIR data). The central amplitude, background am- 
plitude, scale radius, and center position (with respect 
to the center of the SZ signal in SPT data) are free pa- 
rameters in the fit. Priors are placed on 9s, the scale 
radius, and A9, the SZ-NIR positional offset, based on 
the NIR galaxy distributions measured in SPT-selected 
clusters known to be real and typical of the clusters in 
this sample. The number of galaxies above background 
within 1 arcmin^ — which we will call Si/ — is then cal- 
culated from the best-fit profile. The same procedure is 
repeated on fields not expected to contain massive galaxy 
clusters, and the value of Ei' for every SPT candidate is 
compared to the distribution of Ei' values in the blank 
fields. The key statistic is the fraction of blank fields 
that had a Ei' value larger than a given SPT candidate, 
and that value is recorded as Pbiank in Table [3] for every 
high-redshift (z > 0.7) or unconfirmed candidate. This 
technique, including using the blank-field statistic as the 
primary result, is similar to the analysis of WISE data 
in the direction of unconfi rmed Planck Early SZ clus- 
(|20T2|, 



ters in jSayers et al.| ( |2012 p, although that analysis used 
raw galaxy counts within an aperture rather than profile 
fitting. 

The model fitting is performed using a simplex-based 
minimization, with any parameter priors enforced by 
adding a penalty. The positional offset penalty 
is Ax^ — (AO/aAe)'^, where ctas is chosen to be 0.25', 
based on the SZ/BCG offset distribution in Figure [tI and 
the value of r2oo for a typical-mass SPT cluster atnigh 
redshift.'"' A prior is enforced on the scale radius from 
below and above by adding penalties of (0s/^s,max)^ 
and {9s,mm/0sY ^ where 9s.ma.x is chosen to be 0.75'based 
on the 9s distribution in known high-redshift SPT clus- 
ters with NIR data, and 0s,min is chosen to be 0.125'to 
prevent the fitter from latching onto small-scale noise 
peaks. 

For Spitzer /lYiAC and WISE, the fit is performed on 
the 3.6/im and the 3.4/im data, respectively; for NEW- 
FIRM, the fit is performed on the i^s-band data. For 
both Spitzer/IKAC and NEWFIRM, a single magnitude 
threshold is used for every candidate; this threshold is de- 

f'200 is defined as the radius within which the average density 
is 200 times the mean matter density in the universe. 
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termined by maximizing the signal-to-noise on the Ei' es- 
timator on known clusters while staying safely away from 
the magnitude hmit of the shallowest observations. The 
Spitzer/lYlAC data is very uniform, and the 3.6/xm mag- 
nitude threshold chosen is 18.5 (Vega). The NEWFIRM 
/Cg-band data is less uniform, but a magnitude thresh- 
old of 18 (Vega) is safe for all observations. For these 
instruments (IRAC and NEWFIRM), the blank fields on 
which the fit is performed come fro m the Spitzer Deep 
Wide-Field Survey (SDWFS) region ( |Ashby et al |2009| , 
ponds to the Bootes field ot the iNUAU Deep 



which corresponds 



Wide-Field Survey (NDWFS; pannuzi fc Dey"p 99). The 
depth of the SDWFS/NDWFS observations tor both in- 
struments (19.8 in Vega for Spitzer/IRAC 3.6/zm and 
19.5 in Vega for NEWFIRM Ks) is more than sufficient 
for our chosen magnitude thresholds. 

For WISE, in which the non-uniform sky coverage re- 
sults in significant variation in magnitude limits, we per- 
form the blank-field fit on data in the immediate area of 
the cluster (within a ^ 20' radius). Under the assump- 
tion that the WISE magnitude limit does not vary over 
this small an angular scale, we use all detected galaxies 
brighter than 18th magnitude (Vega) in both the cluster 
and blank- field fits. 

3.3. Identifying rBCGs in SPT Clusters 

An rBCG in this work is defined as the brightest galaxy 
among the red-sequence galaxies for each cluster. We 
employ the terminology rBCGs, rather than BCGs, to 
allow for the rare possibility of an even brighter galaxy 
with significant amounts of ongoing star formation, be- 
cause the selection is restricted by galaxy colors along 
the cluster red sequence. We visually inspect pseudo- 
color images built with the appropriate filter combina- 
tions (given the cluster redshift) around the SZ candi- 
date position. We search a region corresponding to the 
projected cluster virial region, defined by 02ooi given the 
mass estimate from the observed cluster SPT significance 
^ and photometric redshift Zp^ot- 

There are 12 clusters out of the 158 with measured 
photometric redshifts that are excluded from the rBCG 
selections. Eight of those are excluded due to contamina- 
tion by a bright star that obscures more than one third 
of the area of the Str SPT positional uncertainty region. 
Another cluster is excluded due to a bleed trail making 
the rBCG selection ambiguous, and three other clusters 
are excluded due to a high density in the galaxy popu- 
lation that, given the delivered image quality, makes it 
impossible to select the rBCG. 

4. RESULTS 

The complete list of 224 SPT cluster candidates with 
SZ detection significance ^ > 4.5 appears in Table [3) 
The table includes SZ cluster candidate positions on the 
sky [RA Dec], SZ detection significance [^], and spec- 
troscopic redshift [zgpcc] when available. For confirmed 
clusters, the table includes photometric redshift and un- 
certainty [zcombifzcomb]: estimated as described in qSTj 
Unconfirmed candidates are a ssigned redshift lower lim- 
its, estimated as described in |3.1.4[ 

We also report a redshift quality flag for each Zphot in 
Table [H For most of the confirmed clusters with reli- 
able photometric redshift measurements, we set Flag = 
1. There is one cluster (SPT-CL J2146-4846) for which 



the three individual photometric redshifts are not statis- 
tically consistent (> 3a outliers) for which we set Flag = 
2. We still report the combined redshift for that cluster 
as in other secure systems. We have 6 cases where we 
only use Swope + Method 1 and 25 cases where we only 
use Spitzer + Method 1 for redshift estimation, marked 
with Flag = 3. There are 2 cases (SPT-CL J0556-5403 
and SPT-CL J0430-6251) where we quote only a Method 
1 redshift even for MOSAIC or IMACS data, marked 
with Flag = 4. In the coadded optical images for SPT- 
CL J0556-5403, we identify an overdensity of faint red 
galaxies at the location of the SPT candidate. This op- 
tical data is too shallow, however, to allow for secure 
redshift estimation, but we are able to measure a red- 
shift by combining this data with NEWFIRM imaging. 
This cluster is the only candidate where we rely on pho- 
tometric redshift from i-K,. SPT-CL J0430-6251 is in 
a field very crowded with large scale structure, making 
redshift estimation difficult. 

In the Appendix, we discuss certain individually no- 
table candidates — such as associations with known clus- 
ters that appear to be random superpositions and can- 
didates with no optical/NIR confirmation but strong ev- 
idence from the NIR overdensity statistic. 
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Fig. 4. — Redshift histogram of 158 confirmed clusters, in bins 
of Az = 0.1. Note that about 18% of the total sample comes from 
z > 0.8. 



4.1. Redshift Distribution 

The redshift distribution for the 158 confirmed clusters 
is shown in Figure [4j The median redshift is z = 0.57, 
with 28 systems (^18% of the sample) lying at z > 0.8. 
The cluster with the highest photometric redshift is 
SPT-CL J2040-4451 at z = 1.35 ± 0.07 (estimated us- 
ing Spitzer/IRAC data) and the highest-redshift spec- 
troscopically confirmed cluster is SPT-CL J0205 -5829 at 
z = 1.32. (This cluster is discussed in detail in Stalder 
e raL][2012| ) 

I'he high fraction of SPT clusters at z > 0.8 is a con- 
sequence of the redshift independence of the SZ surface 
brightness and the arcminute angular resolution of the 
SPT, which is well-matched to the angular size of high-z 
clusters. X-ray surveys, in contrast, are highly efficient 
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at finding nearby clusters, but the mass limit of an X- 
ray survey will increase with redshift due to cosmological 
dimming. ROSAT cluster surveys lack the sensitivity to 
push to these high redshifts except in the deepest archival 
exposures. XMM-Newton archiv al surveys (e.g., Lloyd- 
Davies et al. 2011 Fassbender et al. 2011) and coorHT" 
nated surveys ot contiguo us regions (e.g., Pacaud et al. 



2007 



Suhada et al. 



2012 ) have sufficient sensitivity to 

detect systems like those found by SPT, but the solid 
angle surveyed is curre ntly smaller. For example, the 



Fassbender et al. (2011) survey for high redshift clusters 
will eventually cover approximately 80 deg^ , whereas the 
mean sky density of the SPT high-redshift and high-mass 
systems is around one every 25 deg^ . Therefore, one 



would have expected the Fassbender et al. (2011 ) XDCP 



survey to have found around three clusters ot compara- 
ble mass to the SPT clusters, which is in fact consistent 
with their findings. The vast majority of the high redshift 
X-ray-selected sample available today is of significantly 
lower mass than SPT selected samples, simply because 
the X-ray surveys do not yet cover adequate solid angle 
to find these rare, high mass systems. 

Clusters samples built from NIR galaxy catalogs have 
an even higher fraction of high-redshift systems than 
SZ-selected samples — f or example, in the IRA C Shallow 
Cluster Survey (ISCS; Eisenhardt et al. 2008) a sample 
of 335 clusters has been identified out to z ~ 2, a third 
of which are at z > 1 . However, the typica l ISCS cluster 
M(7, dBrodwin et all 



mass IS 



10^ 



20"07|, significantly 
)Fr high-redshift 



lower than the minimum mass ot the SP 
sample. As with the X-ray selected samples, the Spitzer 
sample includes some massive clusters, including the re- 
cently discovered IDCS J1426.5-H3508 at z = 1.75, which 



was subsequently also detec ted in the SZ ( Stanford et al 



2012 Brodwin et al. 20121). However, the tipitzer sur 



veys to date do not cover the required solid angle to find 
these massive systems in the numbers being discovered 
by SPT. 

4.2. Purity of the SPT Cluster Candidates 

For a cluster sample to be useful for cosmological pur- 
poses, it is important to know the purity of the sample, 
defined as 

so /n\ 

/pure = — ^ , (o) 



where iVtot is the total number of cluster candidates, 
Aureal is the number of candidates corresponding to real 
clusters, and A^faise is the number of false detections. For 
an SZ-selected cluster sample with reasonably deep and 
complete optical/NIR follow-up, a first-order estimate of 
A'lcai is simply the number of candidates with success- 
fully estimated redshifts. In Figure [5] we show two esti- 
mations of purity for the 720 deg^ sPT-SZ sample; the 
first in blue, assuming that all cluster candidates with 
no redshift measurements are noise fiuctuations, and the 
other in red, taking into account incompleteness of our 
follow-up data. The blue/red shaded regions in the fig- 
ure correspond to the la uncertainties on the purity, es- 
timated from Poisson noise on A^faiso for the blue region 
and as described below for the red region. We also show 
the expected purity, estimated from the total number of 
candidates in the sample presented in this work combined 
with the false detection rate from the simulations used 




Fig. 5. — Cumulative purity estimates derived from the opti- 
cal/NIR foUowup compared to simulated purity predictions (black 
line). The inset plot is zoomed-in to the ^ range between 4.5 and 6.0 
and binned more finely. The purity is calculated from the follow-up 
confirmation rate: 1) (blue) assuming all clusters without a clear 
optical or NIR counterpart are false SZ detections (i.e., 100% op- 
tical completeness) and all optical confirmations are robust (100% 
optical purity); and 2) (red) assuming, as justified in the text, 
97% optical completeness and 96% optical purity but taking into 
account clusters confirmed through other means such as X-ray ob- 
servations, la uncertainties in the purity estimates fro m fol low-up 
are shown with shaded blue or red regions (see Section |4.2[|. 



to test the SZ cluster finder (R12 Figure 1). 

The possibility of real clusters beyond the redshift 
reach of our optical/NIR redshift estimation techniques 
makes the blue line in Figure [5] a low er lim it to the true 
purity of the sample. As discussed in { 3.2 we use single- 
band NIR data to estimate the probability that each un- 
confirmed candidate is a "blank field", i.e., a field with 
typical or lower-than-typical NIR galaxy density. Candi- 
dates with no optical/NIR confirmation but with a low 
blank field probability Pbiank, are potential high-redshift 
systems that merit further follow-up study. These sys- 
tems can also give an indication of how much we under- 
estimate our sample purity when we assume any opti- 
cal/NIR non-confirmation is a spurious SPT detection. 
By definition, a low Pbiank implies some NIR overden- 
sity towards the SPT detection, but perhaps not large 
enough to be an SPT-detectable cluster. We can roughly 
calibrate the Pbiank values to SPT detectability by inves- 
tigating the results of the NIR overdensity estimator on 
solidly confirmed, high-redshift SPT clusters. There are 
19 clusters with spectroscopic redshifts above z = 0.7, 
and the average Spitzer/IRAC Pbiank value for these clus- 
ters is 0.04, while the average WISE Pbiank value is 0.05. 
Only three of these clusters have NEWFIRM data, and 
the average NEWFIRM Pbiank value is 0.07. Only one 
cluster in this high-z spectroscopic sample has an IRAC 
Pbiank > 0.1, while three have WISE Pbiank > 0.1. So 
a rough threshold for SPT-type clusters appears to be 
Pbiank ^ 0.1. We havc nine unconfirmed cluster candi- 
dates that meet this criterion in at least one of the NIR 
catalogs, including five that are at Pbiank < 0.05. If we 
assumed all of the Pbiank < 0-05 clusters were real, it 
would imply that the completeness of the optical/NIR 
redshift estimation was ~ 97%, i.e., we have 163 real 
clusters of which we were able to estimate redshifts for 
158. 
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Conversely, the possibility of false associations of spu- 
rious SZ detections with optical/NIR overdensities would 
act in the other direction. Tests of one of the red- 
sequence methods on blank-field data produced a sig- 
nificant red-sequence detection on approximately 4% of 
fields without SPT detected clusters. Assuming that the 
cross-checks with other methods would remove some of 
these, we can take this as an upper limit to this effect. 

We therefore provide a second estimate of purity from 
the optical/NIR confirmation rate, taking into account 
the possibility of real clusters for which we were un- 
able to successfully estimate a redshift (redshift com- 
pleteness < 100%) and spurious optical/NIR associations 
with SPT noise peaks (redshift purity < 100%). From 
the above arguments, we assume 97% for redshift com- 
pleteness and 96% for redshift purity. For each value of 
SPT significance ^, we use binomial statistics to ask how 
often a sample of a given purity with total candidates 
A^(> would produce the observed number of success- 
ful optical/NIR redshift estimates Nconi{> 0) given the 
redshift completeness and false rate. An extra constraint 
is added to this calculation based on data independent of 
the optical/NIR imaging that confirms many of the SPT 
candidates as real, massive galaxy clusters. Specifically, 
we assume an SPT candidate is a real, massive cluster 
independent of the optical/NIR imaging data (and re- 
move the possibility of that candidate being a false opti- 
cal/NIR confirmation of an SPT noise peak) if: 1) it is 
associated with a ROSAT Bright Source Catalog source; 
2) we have obtained X-ray data in which we can con- 
firm a strong, extended source; or 3) we have obtained 
spectroscopic data and measured a velocity dispersion 
for that system. The red solid line and shaded region in 
Figure [s] show the maximum-likelihood value and 68% 
limits for the true purity of the SPT sample under these 
assumptions. 

The purity measured in this work is in good agreement 
with the model for the SPT sample purity. In particular, 
all clusters with f > 6 have identified optical counter- 
parts with photometric redshift estimates. This is con- 
sistent with the expectation of the model and a demon- 
stration that the SPT selected galaxy cluster sample is 
effectively uncontaminated at ^ > 6. With decreasing 
significance, the number of noise fluctuations in the SPT 
maps increases compared to the number of real clusters 
on the sky, and the purity decreases. The cumulative 
purity of the sample is ^70% above — 4.5 and reaches 
^-^100% above ^ = 5.9. Of course, if one requires optical 
confirmation in addition to the SPT detection, then the 
sample is effectively 100% pure over the full sample at 
C > 4.5. 

The high purity of the SPT selected cluster sample is 

jrevious X-ray clus- 



ter surveys (i.e. 


Vikhlinin et al.||1998| Mantz et al. [2008 


Vikhlinin et al. 


2U09), indicating that these mtracluster 



optical follow-up, provide a reliable way to select clean 
samples of clusters for cosmological analysis. 

4.3. rBCG Offsets m SPT Clusters 

The position of the rBCG in galaxy clusters is a prop- 
erty of interest for both astrophysical and cosmological 
cluster studies, as it is a possible indication of a clus- 
ter's dynamical state. In relaxed clusters, it is expected 
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Fig. 6. — rBCG positions are plotted as offsets from SPT candi- 
date positions for 146 systems with clearly defined rBCGs. The ma- 
genta concentric circles enclose 68% (dash-dot line), 95% (dashed 
line) and 99% (dotted line) of the whole rBCG sample and have 
radii of 38.0", 112.1"and 158.6". 



that dynamical friction will tend to drive the most mas- 
sive galaxies to the bottom of the cluster potential well, 
which would coincide with the centroid of the X-ray and 
SZ signatures. On the other hand, in cases of merging 
systems one would expect two different rBCGs, and one 
or both could appear well-separated from the X-ray or 
SZ centroid. Several studies have shown a tight correla- 
tion between the X-ray centroid and the rBCG position 



(Lin&Mohr 
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[[iHaarsm a et al.||2010l [ 

n Fassbcnder et al.ipO] 



jMann & Ebel- 

m 



;liougn J^assbcn aer et ai.| ( |:^Ulip provide ev- 
rat at high redshift the BCG distribution is less 
centrally peaked. Here we examine the rBCG positions 
with respect to the centroid of the SZ signal in the SPT 
cluster sample. 

The position of each rBCG is listed in Table [3j and the 
offsets from the SZ centroids in arcsec are plotted in Fig- 
ure |6] The rings correspond to different fractions of the 
full population of clusters: 68%, 95%, and 99%. These 
rings have radii of 38.0", 112.1"and 158.6", respectively. 
The rBCG population is centrally concentrated with the 
bulk of the SPT selected clusters having rBCGs lying 
within about 1' of the candidate position. 

Given the broad redshift range of the cluster sample, 
the rBCG distribution in cluster coordinates r/r2oo is 
more physically interesting. We use the cluster redshifts 
from this work and the SZ-derived masses from R12 to 
calculate r2oo for each cluster. The red line in Figure[7]is 
the cumulative distribution of the rBCGs as a function of 
r/r2oo- In this distribution, 68% of the rBCGs lies within 
0.17r2oo, 95% within 0.43r2oo and 99% within 0.70r2oo. 

We investigate the importance of the SPT candidate 
positional uncertainty by modeling the expected radial 
distribution in the case where all rBCGs are located ex- 
actly at the cluster center. The la SPT positional un- 
certainty for a cluster with a pressure profile given by a 
spherical /3 model with P = 1 and scale size 6c, detected 
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by SPT at significance ^, is given by 



92 

beam 



(fc^c)2)/C, 



(4) 



where 6'bGam is the beam FWHM , and fc is a factor of or- 
der unity (see Story et al. 2011 for more details). With 
this information, we estimate the expected cumulative 
distribution of the observed rBCG offsets, assuming a 
Gaussian with the appropriate width for each cluster; 
this is equivalent to assuming the underlying rBCG dis- 
tribution is a delta function centered at zero offset with 
respect to the true cluster SZ centroid. Results are shown 
as the blue dotted curve in Figure [7) It is clear that the 
observed distribution of rBCG offeets is broader than 
that expected if all rBCGs were located exactly at the 
cluster center. We conduct a Kolmogorov-Smirnov (KS) 
test to address the similarity of the two distributions. 
The hypothesis that the two distributions are drawn from 
the same parent distribution has a probability of 0.09%, 
suggesting that the observed rBCG offset distribution 
cannot be easily explained by the SPT positional uncer- 
tainties alone. 
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Fig. 7. — Normalized cumulative distribution of rBCG offsets 
from SPT candidate positions as a function of r/r2oo for the SPT 
cluster sample (red line), the expected distribution given SPT po- 
sitional uncertainties if all rBCGs were at exactly the center of 
the true SZ centroid (blue line), and the expected distribution 
given SPT positional uncertainties if the underlying distribution 
of offsets matched those of an X-ray selected cluster sample (green 
line;jLin & Mohr 2004). The KS probability that the observed 
rBCG distribution and the SPT positional error distribution are 
drawn from the same parent distribution is 0.09%, but the ob- 
served distribution is statistically consistent with the distribution 
from the X-ray-selected sample convolved with the SPT positional 
uncertainty distribution. There is no evidence in the rBCG offset 
distribution that SPT-selected clusters are more merger-rich than 
X-ray-selected clusters. 



Because the SPT candidate positional uncertainties IS.9 
are roughly the same, one can expect that our ability to 
measure the underlying rBCG distribution will weaken 
as we push to higher redshift where the cluster virial re- 
gions subtend smaller angles on the sky. We test this by 
dividing the sample into four redshift bins with similar 
numbers of members. The KS tests confirm our expecta- 
tions; using redshift bins of 0.0-0.400, 0.400-0.539, 0.539- 
0.733 and z > 0.733 we find the probability that the ob- 
served and positional error distributions are drawn from 



the same parent distribution is 0.11%, 0.008%, 1.97% 
and 43.4%, respectively. Thus, with the current cluster 
sample, we cannot detect any extent in the rBCG distri- 
bution beyond a redshift z ^ 0.7. If we assume the under- 
lying rBCG offset distribution is Gaussian, the KS test 
shows a maximum probability of 5.3% for a Gaussian dis- 
tributed width of 0.074r2oo with the probability of con- 
sistency dropping below 0.1% for a > 0.08r2oo- There- 
fore, while the Gaussian is not a particularly good fit, 
the measured distribution strongly favors a < 0.08r2oo- 

We test whether our SZ-selected cluster sample ex- 
hibits similar rBCG offsets to those seen in previous X- 
ray studies. To do this, we adopt the previously pub- 
lished BCG offset distribution from the X-ray studies as 
the underlying BCG offset distribution for our sample 
and then convolve this distribution with the SPT can- 
didate positional uncertainties. If rBCGs in SZ-selected 
clusters are no different from those in these previously 
studied samples, then we would expect the KS test prob- 
ability of consistency to be high. We exp lore two sam- 



ples: X-ray model 1 (Lin fc Mohr 2004 green lin e in 



Figurej7| and X-ray model 2 ( |Mann fe Ebeling|2012 



ma- 
genta fine in Figure [7| . The probability ot consistency 
between the SPT sample and X-ray model 1 is 41%, and 
the probability of consistency between the SPT sample 
and X-ray model 2 is 0.37%. 

It appears likely that the differences between the two 
previously published X-ray samples can be explained in 
terms of differences in the BCG selection. The measured 
rBCG o ffset distribution pr esented in this work agrees 
with th elLin fc Mohr| ( |2004[ ) sample, in which the BCGs 
were defined as the brightest if-band galaxy projected 
within the virial radius ^200 with spectroscopic redshift 
consistent with the cluster redshift. This BCG selection 
is very similar to the SPT rBCG selection, with the main 
difference being that we do not have spectroscopic red- 
shifts for all rBCG candidates in the SPT sample. The 
agreement between the SZ- and X-ray-selected samples 
in this case suggests that there are no strong differences 
between the merger fractions in these two cluster sam- 
ples. 



The |Mann &: Ebeling| ( |2012[ ) BCG sample, in contrast, 
was assembled using bluer optical bands, which are more 
sensitive to the star formation history. In addition, in 
cases where a second concentration of galaxies was found 
within the projected virial region, the central galaxy of 
the galaxy concentration coincident with the X-ray emis- 
sion peak was chosen as the BCG, regardless of whether 
it was brighter or not (Mann, private communication). 
This selection criteria would make it difficult to identify 
significantly offset BCGs, which would be more likely to 
be present in merging systems. 

The rBCG offset distribution measured in the SPT SZ- 
selected sample of clusters does not provide any com- 
pelling evidence that SZ-selected clusters differ in their 
merger rate as compared to X-ray-selected clusters. It 
will be possible to test this more precisely once we have 
the more accurate X-ray cluster centers with consistent 
rBCG selections. Currently, the X-ray properties of only 
15 SPT SZ-selected clusters have been published (lAn- 



dersson et al. 
100 additiona 



2011 Benson et al. 20111; however, over 



SP'l selected clusters have been approved 



for observation in on-going programs with Chandra and 
XMM-Newton. With those data in hand we will be able 
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to measure the rBCG offset distribution over the full red- 
shift range of SPT clusters, allowing us to probe for evo- 
lution in the merger rates with redshift. 

5. CONCLUSIONS 

The SPT-SZ survey has produced an approximately 
mass-limited, redshift-independent sample of clusters. 
Approximately 80% of these clusters are newly discov- 
ered systems; the SPT survey has significantly increased 
the number of clusters discovered through the SZ effect 
and the number of massive clusters detected at high red- 
shift. In this paper, we present optical/NIR properties 
of 224 galaxy cluster candidates selected from 720 deg^ 
of the SPT survey that was completed in 2008 and 2009. 
The results presented here constitute the subset of the 
survey in which the optical/NIR follow-up is essentially 
complete. 

With a dedicated pointed follow-up campaign using 
ground- and space-based optical and NIR telescopes, we 
confirm 158 out of 224 SPT cluster candidates and mea- 
sure their photometric redshifts. We show that 18% of 
the optically confirmed sample lies at z > 0.8, the median 
redshift is z = 0.57, and the highest redshift cluster is at 
z = 1.35 ±0.07. We have undertaken a cross-comparison 
among three different cluster redshift estimators to maxi- 
mize the precision in the presented photometric redshifts. 
For each cluster, we combine the redshift estimates from 
the three methods, accounting for the covariance among 
the methods. Using 57 clusters with spectroscopic red- 
shifts, we calibrate the photometric redshifts and uncer- 
tainties and demonstrate that our combining procedure 
provides a characteristic final cluster redshift accuracy of 
Az/(l-hz) = 0.017. 

For the 66 candidates without photometric redshift 
measurements, we calculate lower redshift limits. These 
limits are set by the depth of the optical/NIR imaging 
and the band combinations used. For nine of these can- 
didates there is evidence from NIR data that the cluster 
is a high redshift system, and that we simply need deeper 
NIR data to measure a photometric redshift. 

Under the assumption that all 66 candidates with- 
out photometric redshift measurements are noise fluctua- 
tions, we estimate the purity of the SPT selected cluster 
sample as a function of the SPT detection significance 
^. Results are in good agreement with expectations for 
sample purity, with no single unconfirmed system above 
^ = 6, > 95% purity above ^ = 5, and ^ 70% purity 
for ^ > 4.5. By requiring an optical/NIR counterpart for 
each SPT candidate, the purity in the final cluster sam- 
ple approaches 100% over the full ^ > 4.5 sample. The 
purity of the SPT cluster sample simplifies its cosmolog- 
ical interpretation. 

Next, we examine the measured rBCG offset from the 
SZ candidate positions to explore whether SZ-selected 
clusters exhibit similar levels of ongoing merging as X- 
ray selected samples. We show that the characteristic 
offset between the rBCG and the candidate position is 
~0.5'. We examine the radial distribution of rBCG off- 
sets as a function of scaled cluster radius r jr^m and show 
that a model where we include scatter due to SPT po- 
sitional uncertainties assuming all BCGs are at cluster 
centers has only a 0.09% chance of consistency with the 
observed distribution. That is, the observed distribution 
is broader than would be expected from SPT positional 



uncertainties alone. If we assume the rBCG offset dis- 
tribution is Gaussian, the observations rule out a Gaus- 
sian width of (7 > 0.08r2oo, however, even with smaller 
width a Gaussian distribution is only marginally consis- 
tent with the data. When comparing the SPT rBCG 
distribution with a X-ray selected cluster sample with 
a similar rBCG selection criteria (Lin & Mohr 2004), 
the SPT and X-ray selected rBCG distributions are sim- 
ilar, suggesting that their merger rates are also similar. 
Comparisons to other X-ray selected samples are com- 
plicated by differences in rBCG selection criteria. For 
example, comparing to Mann & Ebeling (20121, which 
selects BCGs using bluer optical bands, we find a sig- 
nificantly less consistent rBCG distribution compared to 
SPT. We conclude that SZ and X-ray selected cluster 
samples show consistent rBCG distributions, and note 
that BCG selection criteria can have a significant effect 
in such comparisons. 

With the full 2500 deg^ SZ survey completed in 2011, 
we are now working to complete the confirmations and 
redshift measurements of the full cluster candidate sam- 
ple. Scaling from this 720 deg^ sample with effectively 
complete optical follow-up, we estimate that the full sur- 
vey will produce ~500 confirmed clusters, with approxi- 
mately 100 of them at z > 0.8. This sample of clusters 
will enable an important next step in cluster cosmologi- 
cal studies as well as the first detailed glimpse of the high 
redshift tail of young, massive clusters. 
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APPENDIX 



NOTABLE CLUSTERS 

SPT-CL J0337-6207 — This candidate is optically unconfirmed but has small NIR blank-fi eld p robability in at least 
one data set (Pbiank = 0.5% in WISE data—also, Pbiank = 18.2% in NEWFIRM data— see p^for details). 

SPT-CL JO428-6O49 — This candidate is optically unconfirmed but has a high SPT signifi canc e — 5.1) and small 
NIR blank-field probability in at least one data set (Pbiank = 0.5% in WISE data, see ^3.2 for details). Though 
^biank — 69.0% in NEWFIRM data, there is clear visual evidence of a NIR overdensity in ttie NEWFIRM images, 
but at ~ 40" from the SPT position. Such a large offset is heavily disfavored by the fitting procedure, such that the 
model that minimizes the overall for the NEWFIRM data is effectively a blank field. The position of the WISE 
overdensity is consistent with the SPT position. 

SPT-CL JO458-574I — This cluster is listed as optically unconfirmed, but it is also listed in Table 2 of R12 as coincident 
with the low-redshift cluster AGO 3298 (at a separation of 77"). We see a clear red-sequence overdensity in our 
Magellan/IMACS data at z ~ 0.2, centered on the Abell cluster position. The best-fit SZ core radius for this 
candidate is 2.5', which implies an SPT positional uncertainty of ~ 0.5', in which case a 77" offset is only a ~ 2cr 
outlier. However, visual inspection of a lightly filtered SPT map shows that the SZ signal is coming from two distinct 
components, one of which corresponds to the Abell cluster position, and neither of which would have been significant 
enough to be included in the R12 catalog on its own. For this reason, we leave the 9c — 2.5' candidate, which blends 
the SZ signal from the two individual components, as unconfirmed. 



This candidate is opticall y un confirmed but has small NIR blank-field probability in at least 
7.5% in WISE data— see P^for details). 



SPT-CL J2002-5335- 
one data set (Pbiank = 

SPT-CL J2032-5627— This cluster is listed in Table 2 of R12 as coincident with the z = 0.06 cluster CIG 2028.3- 
5637 / AGO 3685 (at a separation of 115" from the literature Abell cluster position) and as coincident with the 
z = 0.14 cluster RXC J2032. 1-5627 (at a separation of 87" from the reported REFLEX cluster position). However, 
from our Magellan/IMACS imaging data, we estimate a red-sequence redshift of z = 0.31 ± 0.02, and, using the red- 
sequence measurements as a criterion for cluster member selection, we have obtained spectra for 32 cluster members 
using GISMO and ha ve measured a rob ust spectroscopic redshift of Zgpcc = 0.2840. Examination of the REFLEX 
spectroscopic catalog ( Guzzo et al.||2009 l reveals that their spectroscopic observations yielded five galaxies near their 
reported redshift of ^spec = 0.1380 but also six galaxies within 2% of the value we derive from our GISMO o bservations 
(zsnfir = 0^840). The value of Zspsr = 0-0608 for AGO 3685 is from only one galaxy (and, while reported in Struble & 
is originally from Fetisova 1981). We conclude that there are two clear optical overdensities at ditterent 



Rood 



1999 



redshnts along the line of sight to this SZ/A-ray system, and that the literature redshift of z = 0.0 608 for AGO 3685 



is pro bably incorrect. Because of the redshift dependence of the SPT selection function (see, e.g., Vanderlinde et al. 
20101, it is likely that the bulk of the SZ signal is coming from the higher-redshift cluster. We have obtained XMM- 
iN ewton data on this system, and the X-ray and SZ signals have very similar morphology, indicating that the X-ray 
signal is also predominantly associated with the higher-redshift system. This makes it likely that the z = 0.2840 
system is a massive cluster and that the z ~ 0.1380 system is a low- mass interloper, possibly the cluster originally 
identified as GIG 2028.3-5637 / AGO 3685. SP T-GL J2032-5627 is discussed further— including weak lensing data 
from Magellan/Megacam- 



High et al. (2012). 



This candidate is opticall y un confirmed but has small NIR blank-field probability in at least 
0.1% in WISE data— see P^for details). 



SPT-CL J2035-5614- 
one data set (Pbiank ■ 

SPT-CL J2039-5723 — This candidate is optically unconfirmed but has a small SZ core radius (0.5') and small NIR 
blank-field probability in at least one data set (Pbiank = 1-2% in WISE data and 8.7% in NEWFIRM data— see p!2 
for details). 

SPT-CL J2121-5546 — This candidate is optically unconfirmed but has small NIR blank- field probability in at least 
one data set (Pbiank = 0.9% in WISE data-also, Pbiank = 11-5% in NEWFIRM data— see p^for details). 



This candidate is opticall y un confirmed but has small NIR blank-field probability in at least 
5.2% in WISE data— see sjs^for details). 



SPT-CL J2136-5535- 
one data set (Pbiank - 

SPT-CL J2152-4629 — This candidate is optically unconfirmed but has a high SPT significance (^ = 5.6), a small SZ 
core radius (0.25'), and small NIR blank-field probability in at least one data set (Pbiank = 8.0% in WISE data; also 



Pi 



bla 



10.6% in Spitzer/mkC data and 20.0% in NEWFIRM data— see J 3.2 for details). This is the only candidate 



with Spitzer/mAC Pbiank < 20% for which we were not able to estimate a redshift 



SPT-CL J 23 43- 5 5 56- 
one data set (Pbiank - 



This candidate is optically unconfirmed but has small N IR b lank-field probability in at least 
5.6% in WISE data and 21.0% in NEWFIRM data— see p;2|for details). 
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